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We present an investigation of coherent backscattering of light that is multiple scattered by a 
photonic crystal by using a broad-band technique. The results significantly extend on previous 
backscattering measurements on photonic crystals by simultaneously accessing a large frequency 
and angular range. Backscatter cones around the stop gap are successfully modelled with diffusion 
theory for a random medium. Strong variations of the apparent mean free path and the cone 
enhancement are observed around the stop band. The variations of the mean free path are described 
by a semi-empirical three-gap model including band structure effects on the internal reflection and 
penetration depth. A good match between theory and experiment is obtained without the need of 
additional contributions of group velocity or density of states. We argue that the cone enhancement 
reveals additional information on directional transport properties that are otherwise averaged out 
in diffuse multiple scattering. 



I. INTRODUCTION 

Photonic crystals are ordered nanostructures with a 
periodicity on the scale of the optical wavelength. The 
perfect periodic ordering gives rise to a photonic band 
structure and properties such as slow light, negative re- 
fraction, and photonic band gapsi Photonic crystals are 
of importance for their potential in controlling the trans- 
port and emission of light at the nanoscale^ Apart 
from the properties of perfect photonic crystals, such 
as a 3D photonic bandgap, the possibility of combining 
band structure with controlled random multiple scatter- 
ing has raised considerable interest. Macroscopic pho- 
tonic crystals can be assembled bottom-up using colloidal 
suspensions, or top-down using lithographic techniques. 
All crystals have some degree of intrinsic disorder due 
to polydispersity of the building blocks or imperfections 
in fabrication^— In order to realize random multiple 
scattering without relying on intrinsic disorder, another 
route has been explored to intentionally design disor- 
der in three-dimensional photonic crystals of very high 
quality*^ The role of disorder has recently been investi- 
gated in photonic crystal waveguides i^— 

The presence of disorder results in diffusive scattering 
of light and a concomitant disappearance of the photonic 
band structure. However, theoretical and experimental 
studies indicate that an intermediate regime exists, where 
photonic band structure is preserved locally but medium- 
range disorder disrupts ballistic transport£ i 1Q i 13 i 14 i 17 i 18 . 
This intermediate scattering regime is reached when the 
scattering mean free path £ exceeds the Bragg length 19 
Lb but is smaller than the sample size L: Lb < £ < L. 
As was pointed out by Sajeev John, the modified den- 
sity of states associated with the photonic bands near a 
3D photonic band gap may result in large corrections to 
the diffuse transport, as is known for completely random 



systems £ Recently it has been argued that the scattering 
mean free path in opals with designed disorder yield con- 
tributions from both the local density of states and group 
velocity*^ Ultimately, the photonic crystal could facili- 
tate the complete breakdown of light diffusion known as 
Anderson localization. Indeed, indications of Anderson 
localization have been observed in recent studies on one- 
dimensional photonic crystal waveguide a 15 i 16 and nonlin- 
ear lattices 2 ^. Despite recent experimental efforts^, this 
regime has not yet been observed in a three dimensional 
photonic crystal. 

Effects of multiple light scattering in photonic crys- 
tals have been studied before using a variety of reflec- 
tion and transmission measurements, including coherent 
backscattering of light [Fig. Ufa)]. Coherent backscatter- 
ing (CBS) is a multiple scattering phenomenon caused 
by the constructive interference of time-reversed light 
paths in the medium in a cone around the backscatter- 
ing direction. Coherent backscattering has been used to 
investigate multiple scattering of photonic crystals at a 
few discrete wavelengths A^* 2 ^ Recently the first broad- 
band investigations of opal photonic crystals around the 
stop band have been reported. 22 Modifications of CBS 
cones around the photonic stop bands have been observed 
which were attributed to two contributions related to the 
modification of the diffuse source by Bragg reflection, and 
to variations of the internal reflection conditions at the 
escape angle of the diffuse light, as shown schematically 
in Fig.QJb). Both effects can significantly affect the path 
length distribution for multiple scattering. 

Here, we present experiments demonstrating the effect 
of the photonic band structure on diffuse light scattering 
over a wide spectral and angular range in polystyrene 
opals. We employ a broad-band coherent backscattering 
technique using a supercontinuum white-light source^ 
The technique has been used successfully to investigate 
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FIG. 1: (color online) (a) Experimental setup for coherent 
backscattering spectroscopy, (b) Photonic crystal effects in- 
fluencing the apparent mean free path: 1. modification of 
source of diffuse light by Bragg diffraction, 2. modification of 
angle-dependent escape probability due to internal reflection. 



strongly photonic random media, i.e. Ti02 powder and 
porous GaP, and to characterize resonant light trapping 
by semiconductor nanowires^ Compared to earlier stud- 
ies on photonic crystal o 6 i 7 i 21 i 22 , we achieve a more de- 
tailed characterization of CBS as a function of both an- 
gle of incidence and wavelength. This investigation over 
an extensive range of angles and frequencies allows us to 
compare the results with semi-empirical model calcula- 
tions to verify earlier predictions. 6 



II. EXPERIMENTAL METHODS 

Coherent backscattering spectra were obtained over a 
wide spectral range in the visible and near-infrared us- 
ing a supercontinuum whit e-lig ht backscattering setup 
described previously in Ref. [23|. The angular resolution 
was improved to 0.8 mrad full width at half maximum, 
in order to resolve the narrow CBS cones of photonic 
crystals. 6 Small-angle coherent backscattering measure- 
ments were obtained in a beamsplitter configuration as 
shown in Fig. QJa). Circular polarization channels were 
selected using achromatic polarizers and a quarter wave 
plate directly in front of the sample, with the purpose 
of reducing the contribution of single scattering in the 
CBS enhancement factor. The opals used in this study 
were obtained by sedimentation of polystyrene spheres 
(Duke Scientific, polydispersity ~2%) in 0.3 mm thick, 
3 mm wide flat glass capillaries (Vitro Dynamics). We 
used two opals composed of spheres of radii 130 nm and 
180 nm, which were grown from a suspension in water us- 
ing a centrifuge i 25 i 26 After crystallization the water was 
slowly evaporated from the capillaries to produce poly- 
crystalline opals of several millimeters in length. The 
polystyrene opal is an /cc-crystal with the [111] crystal- 
lographic direction (corresponding to TL) oriented per- 
pendicular to its surface, as was demonstrated by small- 
angle x-ray scattering.— The corresponding Wigner-Seitz 
cell is shown in Fig. Efa), while Fig. [2fb) shows the cal- 
culated band structure and density of states. In our ex- 
periments we investigate the behavior around the stop- 




FIG. 2: (a) Wigner-Seitz cell of an fcc-lattice with its irre- 
ducible part, showing symmetry points T, L, X, U, K, W. (b) 
Calculated band structure for a polystyrene opal, with total 
density of states (DOS) of photonic crystal (line) compared to 
effective medium (dash). Shaded area indicates the reduced 
frequency range accessed in our experiments, with horizontal 
lines indicating individual ranges of the 130-nm (dash) and 
180-nm (dash-dot) opals 



gap at L, which is located at (j) = 0° angle of incidence 
and has a bandwidth of 6% as predicted by band struc- 
ture calculations and measured in many reports on opals 
elsewhere i 12 i 25 i 31 i 5Q i 51 By tilting the sample, we gain ac- 
cess to directions away from the L-gap corresponding to 
LK, LU, and LW in the band diagram. In the experi- 
ments we average over many crystalline regions with ran- 
dom orientations in the (111) plane, and the crystal is 
spun around its TL-axis to average individual speckle. 
Therefore all directions in the (111) plane will contribute 
to our experimental results. Since the LK and LU di- 
rections are very similar, we focus on the LU and LW 
bands in the interpretation of our results. Also shown 
in Fig. [2fb) is the photonic density of states (DOS) 
(line), compared to that of the homogeneous effective 
medium (dashed line). Clearly, density of states varia- 
tions amount to only several percent for opal photonic 
crystals ,2L2£ 



III. RESULTS 

Specular reflectivity measurements were collected for 
different angles of incidence (j). Resulting spectra are 
shown in Fig. [3]^a) and (c) for the r = 130 nm and 
r = 180 nm opals, respectively. For both opals a clear dis- 
persion is observed of the stop band, which is centered at 
16800 cm -1 for the r = 130 nm opal and at 12200 cm -1 
for the r = 180 nm opal at near-perpendicular incidence. 
The frequencies of the half-maximum band edges are 
shown in Figs.[3fb) and (d) together with band structure 
calculations, where the lattice constant was adjusted to 
match the experimental stop band position at (j) = 0°. 
This results in the dimensionless scaled units a/ A as 
shown in the right axes of Fig. [3fb) and (d), and values 
of the lattice constant of 362 nm and 510 nm, respec- 



3 




1.2 1.4 1.6 1.8 2.0 20 40 60 
Frequency (x10 4 cm" 1 ) Angle $ (°) 




1.2 1.4 1.6 1.8 20 40 60 



Frequency (x10 4 cm" 1 ) Angle § (°) 

FIG. 3: Color online, (a) Experimental reflectivity spectra 
of a r = 130 nm opal for angles of incidence from = 4° 
to 50°. Subsequent curves scaled by a decade for visibil- 
ity, (b) Dispersion relations for r = 130 nm opal, (circles) 
half-maximum values of the reflectivity peaks of (a), (lines) 
calculations for fee crystal along LU (dashed lines) and LW 
(full lines) (a = 362 nm). (c,d) Same for r = 180 nm opal 
(calculation, a — 510 nm). 



tively. These values are in good accord with expected 
values based on the sphere radii using a ~ 2y/2r. The 
experimental dispersion of the stop band matches well 
the expected behavior for an fee opal. The full width at 
half maximum of the reflection peaks of Auj/uj is in ex- 
cellent agreement with the photonic interaction strength 
S = 0.06 defined as the relative frequency width of the 
stop gap at the L-point predicted by band structure cal- 
culations. In addition, higher order bands approach the 
stop band at angles <p around 60°, resulting in a more 
complex behavior of the reflected intensity^— , as is ob- 
served for the r = 180 nm sample. Given the spectral 
range covered by our experimental setup, a set of sam- 
ples is needed to span nearly a factor three in frequency 
(a/A ~ 0.35 — 1.05) that ranges from well below the first 
order L-gap to the range of second order gaps that are 
the precursors of the inverse opal band gaps. 

Coherent backscattering spectra were measured for the 
r = 130 nm opal for various angles of incidence </>. A typ- 
ical CBS spectrum taken at (j) = 8° is shown in Fig. Ufa). 
The CBS spectrum shows a pronounced structure around 
16800 cm -1 , corresponding to the center frequency of 
the stop gap. The arrows indicate frequencies below, in- 
side, and above the stop gap, for which cross-sections are 
shown in FigBJb). We note that the height of the CBS 
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FIG. 4: Color online, (a) Experimental CBS spectrum of 
a r = 130 nm diameter opal at 4> = 8° and (b) Cross sec- 
tions at frequencies of 15500 cm -1 (circles, red), 17000 cm -1 
(dots, black), and 18000 cm -1 (diamonds, blue), with (lines) 
fits using diffuse scattering model. (c,d) Values of apparent 
mean free path iw (dots) (c) and CBS enhancement E (d) ob- 
tained from fits to (a). Lines denote cj -0-96 dependence (dash- 
dot, green), obtained from experimental data at (j) — 50° (open 
triangles), and calculations using Eq. (pQ) with (line, red) and 
without (dash, blue) internal reflection correction. 



cones is not at all related to the peak of the reflectivity 
in Fig.[3ja). The traces show the typical CBS lineshape, 
which can be fitted to a model of random light scatter- 
ing [lines in FigHfb)]. We employed a finite slab model 
analogous to Ref. |34|, where the slab thickness was fixed 
for the full data set. The residuals of all fits are less than 
5%, illustrating a very good match. This agreement in- 
dicates that the diffusion Green functions are sufficient 
to calculate the path length distribution underlying the 
CBS cones, without requiring the addition of photonic 
band structure effects. The validity of the random scat- 
tering model is consistent with earlier experiments. 6 The- 
ory has indeed confirmed that isotropic diffuse transport 
is obtained for crystal structures with high symmetry, 
even when scattering between individual modes is highly 
anisotropic^ 

The model fits using the random scattering theory 
yield values for the apparent mean free path tw> which 
relates to the width of the CBS cones according to 
#fwhm — 0.7/ko£w, and the CBS enhancement factor 
E. Spectra of both the apparent mean free path i w 
and the enhancement factor E (obtained from data in 
Fig. Ufa)) are shown in FigHJc,d) (dots). The dashed ver- 
tical lines indicate the high and low-energy edges of the 
photonic stop band obtained from the reflectivity data. 
For comparison we also show the iw spectra obtained at 
(j) = 50° (open triangles). At this angle, the photonic 
stop band has shifted outside of the spectral range, and 
the frequency variation of the mean free path mainly cor- 



responds to a background resulting from structural vari- 
ations in the crystal^ We modeled this variation of the 
mean free path to a power law t oc cj _Q! , yielding an ex- 
ponent a = 0.93 ± 0.04 (dash-dotted line, green). This 
exponent is much less than that of pure Rayleigh scat- 
tering (a = 4), and also less than the Rayleigh-Gans 
dependence (a = 2) found for the transmission of similar 
opals j£ i 13 i 29 i 35 — The low- frequency value of t of around 
20 fim translates to a ratio £/a ~ 55 for our crystal 
(a = 362 nm), in agreement with the long- wavelength 
values observed for similar samples in Ref. l29l . 

For the = 8° spectrum in Fig. IU(c), a pronounced 
trough is observed in £\y at frequencies that correspond 
to the measured L-gap. A similar behavior has been 
observed in experiments^^, and was first predicted in 
Ref. [6| based on a reduction of the penetration depth of 
the diffuse source by Bragg reflection. Qualitatively, a re- 
duced penetration depth implies shorter paths, as paths 
are more likely to meet the sample boundary and exit. 
Hence Bragg reflection causes cone broadening. The re- 
duction of £w extends far outside the stop gap, especially 
toward the lower- frequency side, consistent with the be- 
havior of the reflectivity spectra in Fig. [3fa). This ob- 
servation shows the importance of measuring over a large 
bandwidth in frequency, including a significant range well 
below the stop gap: it is only far below the stop gap that 
one retrieves a transport mean free path that is a true 
measure of the disorder in photonic crystals. The Bragg 
length associated with the L-gap can be estimated from 
the interaction strength according to Lb = 2dm /nS, 
yielding Lb — 4 fim. This corresponds to a ratio 
£/Lb — 3.8, in the L-gap, which fulfills the inequality 
for the intermediate scattering regime (Lb < £ < L). 
The solid and dashed lines in Fig. |4jc) represent model 
fits which will be discussed further below in Sec. IIVDI 

The CBS enhancement factors in Fig. ffl(d) are consis- 
tently low compared to values obtained for strongly scat- 
tering random photonic media measured using the same 
instrument J 3 - This behavior is particularly striking be- 
low the stop band, where photonic crystals are commonly 
expected to behave as an effective mediuni i 41 i 42 The en- 
hancement factor ranges between 1.4 and 1.5, which is 
larger than values reported by Huang et al^ but smaller 
than measured in Ref. @. A reduction of the enhance- 
ment factor below the stop band was also reported in 
Ref. [22|. The ideal enhancement factor of two is expected 
in any multiple scattering sample, whether absorbing and 
of finite size or not, as long as the sample obeys reci- 
procity. However, if scattering components contribute 
that have no distinct time reversed counterpart, such as 
single scattering, the enhancement factor is reduced^ In 
addition to this general trend, we find pronounced vari- 
ations of E around the stop gap. At the band edge, two 
sharp troughs appear in Fig. 0(d), which will be discussed 
further below. 

Results obtained from CBS spectra for the r = 180 nm 
opal at higher reduced frequencies are shown in Fig. [5] 
The spectra for the 180-nm opal yield more information 
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FIG. 5: (a-f) Values of iw (a-c) and E (d-f) obtained from 
fits of CBS spectra for the 180-nm opal. Error bars denote 
95% confidence interval of the fits. Lines: calculations us- 
ing Eq. (pQ) with (line, red) and without (dash, blue) internal 
reflection correction. 



on the behavior above the L-gap, covering values of a/ A 
between 0.55 and 1.03. Values of £w/a are around a fac- 
tor three smaller for this opal than for the r = 130 nm 
sample, which is indicative of an increased disorder. 
Apart from photonic band structure effects, we observe 
a gradual decrease of £ oc uj~ a , which can be modeled 
using a = 0.75 ± 0.08 for the spectrum at <p = 55°. The 
power-law exponent is much lower than that observed in 
transmission measurements, emphasizing the nonuniver- 
sality of the mean free path in opal photonic crystals, as 
opposed to earlier suggestions^ 

Similar to the 130-nm opal, a distinct trough in £w 
is observed in Fig. [5] which shifts with the frequency of 
the stop gap for increasing angle of incidence. From the 
width of the stop gap, we estimate a Bragg length Lb of 
around 5 /im, indicating a ratio £/ Lb of ~ 1 for this more 
strongly disordered opal. We point out the asymmetry 
of both the mean free path and the enhancement factor 
over the stop band in Fig. [51(b). Such asymmetry cannot 
be explained by the group velocity behavior as this is 
symmetrical around the stop gapJ^ 

At angles of incidence above 45° [cf. Fig. [51(c)], a sec- 
ond trough in £w appears around 19000 cm -1 . Compari- 
son with the band structure calculations of Fig. [3] reveals 
that this feature can be attributed to another Bragg re- 
flection at oblique (111) planes, which appears around 
a/A ~ 1 at the U and X points in the Brillouin zone. The 
presence of this band is also clearly observed as a trough 
in the enhancement factor E in Fig. 151(f). The presence 
of this higher-order diffraction peak demonstrates that 
the observed behavior has to be described by the full 3D 
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FIG. 6: Color online, (a) (Left) Color density maps of iw ob- 
tained from fits to CBS cones for the r — 130 nm opal. (Right) 
Calculated iw using Eq. (pQ). (b) Same for the r = 180 nm 
opal. Lines indicate band structure from Fig. [2] (c,d) Calcu- 
lations made using Eq. (pQ). Numbers indicate contributions 
from three bands of Eq. (3]). 



photonic band structure and not just by a simple one- 
dimensional Bragg reflector. 

The extensive information on the diffuse scattering pa- 
rameters obtained from CBS spectra taken at different 
angles and frequencies can be combined into maps of fit- 
ted values for iw and E as shown in Figs. [6] and [71 Fig- 
ures [6ja) and (b) show the experimental data, while Fig- 
ures [6^c) and (d) are results from calculations discussed 
below in Sec. IIV1 The CBS enhancement E is shown in 
Fig. [71 In the following section, we discuss the different 
contributions to the wavelength dependence of iw and 
E in more detail. 



IV. MODEL CALCULATION OF BAND 
STRUCTURE EFFECTS ON 
ANGLE-DEPENDENT CBS CONE WIDTH 

A. General model for CBS-cone width 

In this section we discuss the combined frequency and 
angle-dependence of the width of the CBS cone. Our 
starting point for the description of diffuse light scat- 
tering in the crystal is the theory for random photonic 
media that we extend for the effects of the photonic crys- 
tal band structure According to predictions by Ref. 
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FIG. 7: Color online. Fitted CBS enhancement factor for the 
r = 130 nm opal (a) and the r = 180 nm opal (b). Lines 
indicate band structure from Fig. [2] 



[6|, the cone width is modified by two effects: internal 
reflection 4 ^ and the limited penetration depth of the in- 
cident light when its wavelength matches the Bragg con- 
dition. A larger internal reflection results in longer paths, 
providing a narrowing of the CBS cone. In contrast, re- 
duction of the penetration depth leads to shorter paths 
as the light has a larger probability to exit via the sur- 
face, resulting in a broadening of the cone. It has been 
shown that the resulting correction to the CBS cone can 
be included in the random diffusion model by introducing 
an apparent mean free path iw which is given by£ 



i/i 



w 



(l + r ) 2 C[l + 2( e + r e )C] 
1 + 2t [l + (e + r e )£] : 



with 



R 



l-R 



(1) 



(2) 



Here R is the internal reflection coefficient integrated 
over all angles as put forward in the theory for multiple 
light scattering ] 46 ! 47 The penetration depth of the inci- 
dent light is included in the parameter £ = 1+t/ 'Lb- The 
Bragg length varies between infinity (no Bragg diffrac- 
tion) and 2din/7r£ r , with S = 0.06 the photonic strength 
and dm the lattice constant in the (11 Indirection. We 
use the extrapolation length ratios To = 2/3 and r e = 
To(l + R)/(X ~ R) that were originally derived for ran- 
dom media j 46 i 47 



B. Internal reflection 

We propose to capture the essential optical properties 
of the photonic crystals by the three dominant stop gaps. 
Therefore, the internal-reflection coefficient of the opals 
is modeled as the sum of three Gaussian reflection peaks 



R u(lH) = ^#j(/i;)exp 



(u-Ujbii)) 2 



(3) 
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FIG. 8: Internal reflection coefficient calculated using the 
three-gap model of Eq. ©. 



with angle-dependent peak reflectivities and widths 



Rj(fii) and Aj with jii 



COS( 



Following 



earlier successful descriptions of the angle dependent 
reflectivities^ and escape function 10 , we use the gaps 
along the LU high- symmetry direction to model the re- 
flectivity, and assume a dependence on angle according to 
Rj(fii) = Rj cos (j)+Rj^ g , where the cosine-dependence is 
based on simple Bragg diffraction. The frequencies 001,2,3 
are obtained from the calculated band edges in Fig.[3fd). 

Using the three-gap model described by Eq. (j3j), the 
internal reflection R is obtained by angle- averaging for 
every frequency using 4 ^ — 



R L . 



3G 



2,u 



2C 1>a 



(4) 



with 



tt/2 



R®(a) cos n (a) sin a da. 



(5) 



The resulting values of the internal reflection for an opal 
photonic crystal are shown in Fig. [8j The internal reflec- 
tion due to photonic band structure ranges between zero 
and 14% in the frequency range of interest. We notice 
three peaks in the spectrum of R, which can be attributed 
to the L-gap and the avoided crossing at the U and W 
points. The internal reflection contribution should not be 
confused with effects due to density of states, which also 
are independent of angle. There is a fundamental physi- 
cal difference between these two physical parameters: the 
internal reflection coefficient is the result of gaps, in other 
words, of the absence of states. Conversely, the density 
of states is a sum over states, and is thus the result of 
the presence of states. Moreover, the variations in R 
are much greater than those in the density of states, see 
Fig.Etb) and Refs. [H,[27|. 



Modification of diffuse source by Bragg 
diffraction 



At the center of the stop gap, the following expression 
holds for the Bragg length in a two-band model 4 ^ 



Li 



2dm 
ttS 



(6) 



The Bragg length varies between Eq. © at the Bragg 
condition and infinity outside the stop band. In absence 
of a theoretical description of Lb , we assume here Lb to 
be inversely proportional to the reflectivity around the 
stop gap 



L B (m,oj) = 



2^ 



11 



1 



ttS RS(jh)' 



(7) 



where R® (/i^) is the internal reflection coefficient defined 
in Eq. (|3j) and S = 0.06 is the photonic strength obtained 
from band structure calculations and experiments. This 
definition of Lb is chosen as it provides the right value 
in the middle of the gap, the peak in 1/Lb follows the 
reflection coefficient as expected, and outside the reflec- 
tion peaks Lb becomes very large, as expected. Further- 
more, a reduced reflection coefficient yields an increase of 
Lb-, which is qualitatively consistent with experiments by 
Vlasov et al. and Bertone et a h 44 i 49 More accurate mod- 
elling of Lb is beyond the scope of our semi-empirical 
analysis. 

The modification of the diffuse source by Bragg re- 
flection is included in Eq. (pp) through the parameter 
£ = 1 + I /Lb- The parameter £ can be identified as 
the inverse penetration depth of the light in the crys- 
tal due to Bragg reflection^ Below we will calculate the 
effect of photonic stop bands on £ for the opals under 
study. We note that the model presented for internal re- 
flection here is a rigorous framework in multiple scatter- 
ing theory, barring the introduction of a semi-empirical 
description of the band structure in Eq. (j3j) and Eq. ([7]). 
These approximations limit detailed comparison particu- 
larly at the edges of the stop band. Our model is highly 
constrained, as the only adjustable parameters that we 
allow to vary are the ^1,2,3 that are constant or very 
smoothly varying. 



D. Comparison with experiments 

Using the model above, we will estimate the effects 
of the internal reflection and source distributions on the 
CBS cone width. Values of £, R±, and R2 were optimized 
to achieve good agreement with the experimental data 
of Fig. EJ yielding R x = 0.5 + 0.2 cos 0, R 2 = 0.02, and 
Rs = 0.4. These four coefficients together with the over- 
all scaling of the mean free path are the only adjustable 
parameters used to fit the data for all frequencies and 
incidence angles. The small value of R2 is chosen since 
this band does not have a pronounced appearance in our 
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FIG. 9: Color online, (a) CBS cone narrowing factor £w/£ calculated using Eq. (pQ) including effect of internal reflection R 
from Fig. [8] (b) Calculated penetration depth parameter £, and (c) CBS cone narrowing resulting from combined effects of R 
and £. 



experimental data. Figure [9fa) shows the effect of in- 
ternal reflection only on £w- The increase of R results 
in a narrowing of the cone and thus an increase of £w 
according to Eq. (pQ), which does not depend on angle of 
incidence <\>. In contrast, the inverse penetration depth £, 
shown in Fig. [9jb), follows the band structure according 
to Eq. ©. 

The two contributions are combined to yield the ap- 
parent mean free path as shown in Fig. [6fc,d) for the 
two opals. For comparison with the experimental data, 
we included the power law dependence of the background 
mean free path in absence of photonic effects. The calcu- 
lated maps for £w in Fig. [6] agree well with most of the 
experimental observations. The cross sections, shown as 
red lines in Figs. 2] and [5j reveal that the depth of the cor- 
rection on £w closely matches the experimental behav- 
ior. However, the effect of internal reflection, predicted 
as an increase of £w before the stop band at larger an- 
gles, is not consistent with our experimental data. We 
note that effects of internal reflection are generally diffi- 
cult to disentangle from smoothly varying corrections on 
£w of the order of 10%. The predominant modification 
of the cone width in the photonic crystals under study is 
caused by the reduction of the penetration depth, as is 
illustrated by calculations including only this correction 
(dashed lines, blue). 



E. Enhancement factor 

Here, we will briefly discuss the remarkable behavior of 
the CBS-enhancement factor. After elimination of spuri- 
ous background contributions to the intensity, variations 
in E can only be caused by changes in the intensity and 
polarization of single-scattered light. The low enhance- 
ments at small reduced frequencies are particularly puz- 
zling, since photonic crystals in the long wavelength limit 
are widely expected to behave as effective weakly scat- 
tering media i 25 i 41 i 42 A similar reduction of E in the long- 
wavelength limit was observed for thin slabs of nanowires, 



which was explained by the finite thickness L of the layer 
of the order L/£ < 10 and the concomitant suppression 
of higher scattering orders^ For the current samples, in 
the long wavelength regime, L/£ is still larger than 15. 
Therefore we can exclude finite sample thickness as an 
origin for the low enhancement. 

The effect of stop gaps on single scattering may be 
caused by changes in intensity or polarization of the 
single-scattered light. While absorption loss does not 
change the enhancement factor E = 2 for an ideal ex- 
periment where only multiply scattered light is collected, 
absorption might influence E via the ratio between single 
and multiply scattered light. However, absorption would 
also be apparent as a cone rounding. By analyzing the 
cone rounding, we note that the diffuse absorption length 
is equal to or even greater than the sample thickness, 
i.e. 200 /im. This value corresponds to an absorption 
decay length of 6 mm, in agreement with our previous 
measurements. 51 Therefore, we conclude that the role of 
absorption loss is negligible. 

Variations in the total scattering intensity are observed 
around the stop gap in our samples, as shown in Fig. [TUJ 
consistent with the behavior reported in other work^^ 
We see a minimum in the stop gap with two wings of 
high scattering at the band edges. Inside the stop gap, 
the Bragg reflection removes intensity from the diffuse 
scattering. This redistribution does not necessarily influ- 
ence the CBS enhancement factor considering that the 
Bragg peak does not lead to additional single scattering. 
A decrease of the incoherent single scattering background 
in the stop gap may even be expected as single scattering 
from the photonic crystal coherently adds to the specular 
Bragg reflection and is thus suppressed at other angles. 
The overall increase of scattering intensity toward shorter 
wavelengths and larger angles of incidence is consistent 
with an increase of the enhancement factor which can be 
attributed to more efficient randomization of light in the 
medium. 

Another pronounced feature in E are the troughs at 
the band edges observed for the r = 130 nm opal in 
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FIG. 10: Scattering intensity from the 130-nm (a) and 180- 
nm (b) opals, for angles of incidence from 8° to 50° , measured 
at a detection angle of 5° using helicity conserving channel. 
Subsequent spectra are offset by 0.5 for visibility, starting 
from <j) — 50° (no offset). 



Fig. The conditions of these throughs appears cor- 

related with presence of sharp peaks in the scattering 
intensity of Fig. [TUt a), which currently is not explained 
by any models. These effects resemble the characteristic 
behavior of the group velocity as well as recent measure- 
ments of scattering mean free path^ Such a reduction 
of the single scattering length £ s at the band edges is not 
observed in the apparent mean free path £w of Fig.H^c), 
as £w is averaged over many different directions in the 
crystal. We propose that the CBS enhancement factor 
is inherently more sensitive to directional band struc- 
ture effects through its dependence on single-scattering, 
and therefore changes of £ s can result in variations of 
E. The large variations in E observed in our experiment 
call for new theory for diffuse light in (partially) ordered 
nanophotonic materials. 

Some of the variations in the CBS enhancement cannot 
be associated with the behavior of the total scattering in- 
tensity in Fig. [TUJ but may be associated to the use of 
circularly polarized light in our experiments. Markedly, 
there is the drop in E associated with the second stop 
band at (j) = 55°. This opposite behavior to the increase 
of E in the L-gap points to a difference in helicity of sin- 
gle scattering, possibly due to Brewster angle effects. 52 
For the L-gap at </> = 0°, the parallel and perpendicular 
components perceive exactly the same Bragg reflectivities 
and thus no change of the polarization state of circularly 
polarized incident light should occur. Strong polariza- 
tion dependence of light transmitted and reflected from 



opals was identified in several studies j 5Q i 52 i 53 Polariza- 
tion effects have also been found in Ref. [22|, where in 
particular cross-polarized CBS cones were reported. Our 
measurements of the cross-polarized channel only showed 
an angle-independent response; no cross-polarized CBS 
cones were detected. 

V. CONCLUSIONS 

We have investigated multiple scattering of light in 
three-dimensional polystyrene opal photonic crystals us- 
ing white-light coherent backscattering spectroscopy. 
Our experiments significantly extend on earlier measure- 
ments by covering simultaneously a large frequency and 
angular range. The measurements of the backscattering 
cone width at frequencies near the stop bands confirm 
the model by Refs. [6.45] which includes internal reflec- 
tion and the limited penetration depth of the incident 
light. The semi-empirical model matches the experiments 
without the need to include effects of density of states or 
group velocity, such as were proposed in Ref. [13| . In ad- 
dition, we report photonic band structure contributions 
in the CBS enhancement factor. We propose that the 
CBS enhancement factor, being particularly sensitive to 
single scattering, thus contains information on properties 
which are washed out in other aspects of the multiple 
scattering CBS cone. A better understanding of group 
velocity and other band structure effects to the enhance- 
ment factor will require further theoretical studies. 

The presented experimental approach using white- 
light coherent backscattering opens up studies on diffuse 
light transport in novel materials. Future investigations 
may look toward stronger scattering photonic crystals 
such as Ti02 or silicon inverse opals or silicon inverse 
woodpiles 5 - 5 -, for which effects of the density of states 
could be observable. Other examples where the interplay 
of disorder - or reduced symmetry - and band structure 
is expected to result in novel phenomena are photonic 
quasicrystals 5 -^ and optical graphene j 57 i 58 
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